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ABSTRACT 

A composite method for the automatic detection of diagnostic features related to the depolarization sequence (P-QRS complex) of 
the heart, for arrhythmia classification, using single lead ECG is presented. The non-syntactic approach based upon slope and 
amplitude thresholds along with a set of empirical criteria is employed for segmenting QRS complexes from a variety of noisy ECG 
recordings acquired from the MIT/BIH arrhythmia database. The background noise is removed from the non-QRS portions using an 
appropriate filtering method that causes no change in the amplitudes or boundaries of P and T waves. In R-R intervals, the isoelectric 
line is determined by developing a technique based upon the method of least squares approximation. Amplitude threshold bands are 
set up above and below the isoelectric line to detect P and /or T wave peaks. Using a group of decision logic rules, framed on the 
basis of an exhaustive study of normal and arrhythmic ECG signals and detailed consultations with two independent cardiologists, 
the P waves are discriminated from the T waves. In total, 37 useful diagnostic features have been deduced pertaining to the QRS and 
P waves in the time domain. The QRS detection algorithm of the composite method was validated using 32,800 beats of several 
records of the MIT/BIH database for which a detection accuracy of 99.96% was achieved within the tolerance limits recommended 
by the CSE Working Party. The composite method for detection of both P and QRS wave features has been validated using all the 25 
records for lead II of the CSE Dataset-3, before being applied to approximately 8000 beats of the MIT/BIH database. As regards 
noisy signals, only those were analysed that had a baseline wander not exceeding 0.25 Hz. 


Keywords: Arrhythmia, ECG signal, P-QRS complex, MIT/BIH data base, atrio-ventricular conduction ratio, diagnostic feature 


INTRODUCTION one complete ECG cycle or beat demonstrating the 
wave intervals and segments normally observed in the 
Cardiac arrhythmias occur when disturbances are ECG is shown in Fig. 1. 

caused in the normal electrical events related to the 

A normal ECG rhythm is an ordered sequence of 
depolarisation of the myocardial cells [3], ie. the 


sequence of P wave (atrial myocardial depolarisation) 


basic processes of automaticity, conduction and 
triggering mechanisms of the heart [1, 2]. If not well 


diagnosed in time, they represent a serious threat to the ; : foi, 
and QRS wave (ventricular myocardial depolarisation) 


patient. Therefore, there is a need for quick . : 

: sien il : a: generation, at a regular rate of 60-100 beats per minute 

identification of these abnormal electrical activities of : : . Pa 
(bpm). When this sequence is disturbed, it is 


the heart. Both life threatening (e.g. ventricular manifested as an irregularity in the heart rate or ECG 
fibrillation and atrial fibrillation) and not-so-life morphology. As a result, abnormal rhythms or 
threatening arrhythmias (e.g. premature ventricular arrhythmias occur. 


contraction and atrial premature contraction) can be 


detected by using the electrocardiogram (ECG) that is Rhythms analysis to diagnose arrhythmias involves the 


a low-cost, non-invasive and standard diagnostic 
measure employed by cardiologists all over the world. 


A standard ECG record is the best test for diagnosing 
all arrhythmias, whether of ventricular or 
supraventricular origin. An ECG tracing is a series of 
waves that represent the electrical events of the various 
chambers and conduction pathways within the heart. 
The electrical activity during the cardiac cycle is 
characterized by five separate waves of deflections 
designated as P, Q, R, S and T. The configuration of 


accurate detection of P waves and QRS complexes 
with respect to time, with respect to each other and 
with respect to space [3, 4]. Arrhythmia monitoring is 
performed generally using ECG lead II, since both the 
P waves and QRS complexes are observed clearly in 
this lead [1-6]. 


Computer based rhythm analysis is the state-of-the-art 
technology assisting cardiologists and physicians now- 
a-days. A large number of algorithms have been 
developed by researchers to combat the problems 
cardiologists face while performing the tedious task of 
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visual interpretation of Holter ECG records. The work 
carried out in this paper is another attempt at 
automated extraction of diagnostic features from ECG 
signal that would help the cardiologist to perform 
reliable arrhythmia detection. 


A composite methodology has been developed in the 
time domain that extracts 37 useful diagnostic features 
from ECG signal which is an improvement over 
previously reported algorithms [7-23] used for 
analysing the ECG signals of the MIT/BIH standard 
arrhythmia database. 


CHARACTERISTICS OF 
ARRHYTHMIA DATABASE 


MIT/BIH 


The MIT/BIH database consists of roughly 109,000 
beats that have been manually annotated by at least 
two cardiologists working independently. The 
reference annotation files, include beat, rhythm, and 
signal quality annotations. Each signal file contains 
two signals sampled at 360 Hz. 


Modified lead II (MLII) has been provided as one of 
the two channel ECG recordings [24]. Most of the 
arrhythmia detection algorithms reported in literature 
have opted MLII for performing QRS detection and 
rhythm analysis [9, 15, 18, 25]. However, ECG record 
segments, with poor signal quality and corrupted by 
high noise have generally been excluded from rigorous 
analysis due to the deteriorating performance of the 
QRS detection software while processing these records 
to accomplish arrhythmia detection [6, 10, 13, 18, 20, 
24, 26]. 


It is to be noted that real time ECG signals of the 
MIT/BIH database are largely contaminated by 
different frequencies of noise, viz., baseline wander in 
the range of 0.05 to 0.8 Hz, most of the signals having 
a baseline drift well below 0.5 Hz [6,12,26-29]; motion 
artifacts below 7 Hz [12]; powerline interference at 60 
Hz, and multiples of 30 Hz [24]; and other high 
frequency noise, like muscle noise in the range 50 Hz 
to 2 kHz [18,30,31]. 


The most useful diagnostic components of the ECG 
signal, for arrhythmia analysis, the P wave and the 
QRS complex, also lie within these frequency ranges. 
The relative power spectra of QRS complex lies 
approximately within 5 to 25 Hz, the maximum QRS 
energy being around 17 Hz, while P and T waves lie in 
the frequency range below 10 Hz [6, 12, 25, 26, 32]. 
Therefore, while removal of the multiple frequency 
noise, it is very important to decide how much to filter, 
i.e. how to preserve the “true” signal content to be used 
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for diagnostics, while eliminating noise that interferes 
during the process of feature extraction [18,28,33]. 


Another important aspect, related to the annotated 
ECG records of the MIT/BIH database is that the 
database does not provide any measurement data that 
can be utilized to verify the magnitudes and boundaries 
of different features extracted (e.g. QRS-onset, QRS- 
offset, P-onset, P-offset). 


Manual measurements or seeking help from a 
cardiologist are the only ways to verify the extracted 
QRS and P wave features. In this scenario, an attempt 
has been made to develop a method for automated 
feature extraction in the time domain using real-time 
MLII ECGs of the MIT/BIH database. 


However, prior to rigorous testing and application of 
the QRS and P wave detection algorithms using the 
records of the MIT/BIH database; these algorithms 
were validated using the 25 records of the CSE 
Dataset-3 for which the measurement data are provided 
[34]. 


PRESENT 
PROCESSING 


METHODOLGY FOR ECG 


Composite Method of Feature Extraction 


The technique developed performs cycle-by-cycle or 
R-R interval by R-R interval feature extraction from 
the given real-time ECG signal to enable interpretation 
of different arrhythmias. Cycle-by-cycle feature 
extraction helps the cardiologist while diagnosing one 
arrhythmia (e.g. sinus bradycardia, atrial fibrillation) at 
a time, or a combination of two or more arrhythmias 
(e.g. sinus bradycardia with right bundle branch block) 
coexisting in the signal. 


Also, this approach aids the cardiologist in detecting 
multiple arrhythmias in the same sequence as they 
appear in the ECG signal under examination. 


The composite methodology developed for automated 
extraction of both QRS and P wave features for 
arrhythmia classification, comprises of following 
major elements: 


e ~=Extraction of QRS complex from noisy non-QRS 
segments in the ECG signal; extraction of QRS 
complex features such as _ peak (Q,R,S) 
amplitudes, QRS onset and offset etc.; derivation 
of associated features such as R-R interval, heart 
rate, QRS configuration etc 


e = Digital filtering of noisy non-QRS segments 
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e = Detection of iso-electric (IE) line in R-R interval; 
IE line is used as reference base line to detect 
peaks present in non-QRS segments in R-R 
interval 


e Extraction of P wave and associated features from 
non-QRS segments of R-R_ intervals using 
decision logic rules (employed to discriminate 
between P and T waves); detection of features 
such as P wave position and polarity, P wave onset 
and offset, and derivation of associated features 
such as P-R interval and P-P interval, etc. 


The complete procedure implemented, in the present 
study, is outlined in Fig. 2. The QRS detection 
algorithm, developed in the present study, first 
identifies the R peak point in the QRS complex. Using 
the detected R peak point as a reference, the algorithm 
is further extended to detect other component peaks of 
the QRS complex, QRS onset and offset and other 
associated diagnostic features. 


The detection of P waves in the ECG signal involves a 
pre-processing stage that eliminates noise and 
smoothens the non-QRS segment present in R-R 
interval. Thereafter, the IE line is computed for each 
R-R interval to account for abrupt changes in baseline 
of the ECG signal. 


A point of inflexion in smoothened non-QRS segment 
is taken as a P or T peak point if it falls within the 
amplitude threshold band limits defined above and 
below the computed IE line. The P peak is 
discriminated from T peak using a set of decision logic 
tules developed in the present study. The technique 
eventually computes the diagnostic features associated 
with the P wave. Important details of the composite 
method are as follows: 


ECG Segmentation — Extraction of QRS Complex 


QRS complexes are extracted from the noisy ECG 
signals using a method developed, based on the very 
widely employed non-syntactic approach: _ the 
amplitude and slope threshold technique [13,14, 18, 
30, 32, 35-38] that is immune to low-amplitude noise 
corrupting the non-QRS portions of ECG signals, and 
provides successful QRS detection in the presence of 
moderate baseline wander. 


The frequency spectrum of the QRS complexes 
overlaps that of the noise contaminating non-QRS 
segments. Hence, QRS complexes are segmented out 
from the noisy ECGs and analysed separately from the 
remaining non-QRS noisy portions of the ECG signal 
that need to be filtered before further processing to 
detect P waves. 
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The background noise is rejected by QRS detection 

procedure that is based on computing the slope that lies 

within a certain threshold. The slope is computed using 

first difference for an optimized number of seven 

consecutive samples along the rising R wave, and the 
change in its sign is observed. 


At the point of inflexion, sign of the slope changes 
from positive to negative (or vice versa). This point of 
transition is taken to be a possible QRS complex peak, 
however, its final detection is subject to certain 
additional criteria consisting of decision rules 
providing optimum thresholds that characterise various 
QRS wave features. 


These criteria are formulated in the present 
methodology on the basis of standards of ECG 
interpretation [1-3, 26, 35, 37], recommendations of 
the CSE Working Party [39], statistical variation of the 
ECG data processed in the present study, and extensive 
trial-and-error tests conducted over a large quantity of 
ECG data. 


These additional criteria are as follows: 


(a) The general configuration of the QRS complex 
(with upright R wave), in the MLII ECG signals, 
contains Q-R-S-R’ peaks in a sequence. Therefore, 
a detected QRS complex may be a peak complex 
[39, 40], having the following morphologies: 
QRS, QR, RS, R, RSR’, QRSR’, QS (Table 1). 
Annotated beats in the MIT/BIH database with 
single negative waves have been termed in the 
present study as QS waves, 


(b) QRS onset has been used as the reference point for 
determining the base level to perform all 
amplitude measurements pertaining to QRS 
complex [39], 


(c) R’ peak cannot exist without S peak in a QRS 
complex [36, 40], 


(d) ECG signals with acceptable noise levels can be 
subjected to noise thresholding to detect Q, S and 
R’ peaks in presence of noisy peaks of low 
amplitude [37, 40, 41]. After extensive tests 
conducted over a large number of ECG data, an 
optimized amplitude threshold of 0.06 mV was 
determined for differentiating Q, S and R’ peaks 
from noisy ECG peaks (Daskolov et al. [37] have 
used an amplitude threshold of 0.05 mV to 
differentiate the component peaks of QRS 
complex from noisy peaks. The difference 
between the amplitude thresholds arrived at in the 
present study and that of Daskalov et al. may be 
attributed to the fact that their study was carried 
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out using a 12 lead recorded data instead of the 2 

lead MIT/BIH database used in the present study), 

ie.: 

e Any negative peak before the R peak, with an 
amplitude of 0.06 mV or greater with respect 
to (w.r.t.) the QRS onset, is the Q peak 


e Any negative peak after the R peak, with an 
amplitude of 0.06 mV or greater w.r.t. QRS 
offset, is the S peak 


e Any positive peak after the S peak, with an 
amplitude Of 0.06 mV or greater w.r.t. the 
QRS offset, is the R’ peak 


(e) The noise corrupted nature of ECG signal is 
utilized for detecting onset and offset of QRS 
complex. Onset of Q wave (if present) is the onset 
of QRS complex at a point where the slope 
threshold of QRS wave changes or the slope 
becomes nearly zero [35, 37]; onset of R wave is 
the onset of QRS complex if Q wave is absent. 
Similarly, offset of QRS complex is the offset of 
R’ wave (if present) or is offset of S wave when 
R’ wave is absent or is the offset of R wave when 
S wave is absent, 


(f) A slur or a notch may occur at any sample point 
either along the rising or falling limb of the 
component waves (Q, R, S and R’), interfering 
with computation of onsets, offsets and peak 
extremums of these waves. Fig. 3 shows eight 
types of notched and slurred wave configurations, 
where A, C, E and G are slurs, and, B, D, F and H 
are notches that may be encountered in the QRS 
complexes. To deal with this practical situation a 
suitable method is employed to test for the QRS 
boundaries, slurs and notches. Mainly, a check is 
made for the calculation of the respective peak 
points, offsets and onsets. For instance, if a 
deflection (C or D type) occurs in the rising limb 
of the S wave, then for next few data samples 
beyond this deflection, it is checked whether the 
upward slope continues within a threshold or not. 
If it continues, then the deflection encountered in 
the wave is specified a slur /notch, and the search 
is further carried out to detect the offset of S wave 
or peak extremum of R’ wave. If the upward slope 
does not continue within a threshold for a few data 
samples beyond this deflection, then the first 
sample point of the occurrence of this deflection is 
marked as the offset of S wave, and hence the 
offset of QRS complex. In a similar manner, other 
types of slurs and notches may be distinguished 
from the onsets, offsets and peak extremums of Q, 
R, S and R’ waves. This technique was tested on 
slurs and notches of small magnitudes, 
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(g) In more than 90% of the QRS complexes of the 
MIT/BIH database records (having a_ baseline 
wander in the range of 0.05 Hz to 0.25 Hz), 
analysed in the present study to calculate QRS 
duration, it was found that the amplitudes of QRS 
onset and offset differed within a range of 0.1 mV. 
Therefore, this range was utilized to detect 
different features of the QRS complex such as 
QRS offset, S peak and R’ peak in presence of 
notches and slurs. Results of the computations of 
the QRS duration measured from its onset to 
offset, after implementing the above condition, 
were found to be within the prescribed tolerance 
limits [39], 


(h) In case a T wave occurs with a magnitude 
comparable to that of the preceding QRS complex, 
the criteria used by Pan and Tompkins [13] is 
employed to discriminate between the T wave and 
QRS complex. 


The step-by-step procedure developed for detection of 
QRS complex and associated features is illustrated in 
Fig. 4. 


Digital Filtering of non-QRS Segments 


The non-QRS portions of the MLII ECG signals are 
corrupted by noise. This noise was removed to identify 
P and T peaks, and subsequently extract only P waves 
using the technique of amplitude and slope threshold 
criteria. For computation of slope, the simplest 
technique of first difference has been used that 
performs very well for noise free signals [32]. 


In the present study, a 6""-order low pass Butterworth 
digital filter having a cut-off frequency of 10 Hz is 
used to remove the noise in the non-QRS segments of 
the acquired real-time ECGs. Its amplitude response 
with cut off frequency of 10 Hz is shown in Fig. 5. 


This 6'-order IIR (infinite impulse response) digital 
filter is a cascade combination of three 2"’-order filter 
sections, each 2" order section generating a separate 
set of five digital filter coefficients that are used to 
produce a noise free signal. The filter design is 
accomplished by the algorithm developed by Stearns 
and David [42]. The primary advantages of this filter 
are its sharp frequency cut-off and stability [6, 42]. 


To select an optimum cut-off frequency for filtering 
noise in the non-QRS segments by the low pass 
Butterworth filter, Fast Fourier Transform (FFT) was 
computed for a large set of these segments acquired 
from different records of the MIT/BIH database. 
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From these FFT responses, it was found that the 
frequency content of P and T waves lies within 10 Hz. 
Hence, 10 Hz was selected as the optimum low-pass 
filter frequency cut-off. It was also observed that this 
low pass filter scheme provides noise free and 
smoothened signal with no substantial amplitude 
attenuation of P and T peaks and no change in the R-R 
and P-P intervals. 


Detection of IE Line in R-R Intervals 


The imaginary JE line in an R-R interval is determined 
between two consecutive QRS complexes of the 
filtered ECG. It is used as a reference base line to 
detect P and T peaks. However, only P wave 
parameters have been extracted from the non-QRS 
interwave segments of R-R intervals. 


Based on the recommendations to determine a 
reference baseline [39], a technique has _ been 
developed to determine an IE line in an R-R interval. 
However, it is to be noted that in the arrhythmic ECG 
signals of the MIT/BIH database, the P wave may fuse 
with QRS onset or may not be present at all (e.g. 
Premature Ventricular Complex- PVCs, ventricular 
tachycardia etc.) in some portions of the signal under 
analysis. 


In such cases, the P-R segment criteria for determining 
the IE line cannot be implemented. Therefore, the 
technique determines an IE line that passes through 
samples within a certain amplitude threshold using the 
onset of QRS complex as a reference point. These 
samples lie in the window having a maximum length 
of 66% of R-R interval starting from the onset of QRS 
complex. 


Macfarlene et al. [43] used linear interpolation to 
determine the baseline in ECG signals. In this study 
too, a straight line is fitted using the technique of least 
squares approximation [42] by scanning the samples 
located in the window. The onset and offset of this 
imaginary IE line are then, respectively, compared 
with the detected offset of previous QRS complex and 
onset of current QRS complex. 


If the comparison satisfies a stipulated threshold (that 
was arrived at after extensive trial and error tests on a 
large number of ECG data with moderate baseline 
wander lying within 0.05 to 0.25 Hz), then the fitted 
imaginary IE line is taken to be finally detected, 
otherwise the window size for the R-R interval is 
modified (in terms of both the length and height of the 
window). The window size is modified (i.e., the 
window size contracts or expands sequentially) a 
maximum of four times in order to detect a horizontal 
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/slant imaginary IE line. The step-by-step procedure 

for detection of JE line between two consecutive QRS 
complexes is illustrated in Fig. 6. 


P Wave Extraction from ECG Signal 


In most of the arrhythmia monitoring systems the next 
step after valid QRS complex recognition is P wave 
detection followed by extraction of P wave features. 
Extraction of low amplitude P waves in a noise 
corrupted ECG has posed problems in_ several 
arrhythmia detectors. Therefore, most arrhythmia 
monitors used in hospitals and clinics depend on 
features related to QRS complex only, like R-R 
interval, heart rate, QRS duration etc. [2, 3, 6, 31, 35, 
44]. 


However, in the present study, P wave detection (based 
on the amplitude and slope threshold criteria) has been 
successfully performed, whatever is its magnitude, 
polarity or location in an R-R interval, since, non-QRS 
segments are filtered and smoothened prior to P wave 
extraction. The methodology for detection of P waves 
and associated features is illustrated in Fig. 7. 


Amplitude threshold band limits above and below the 
computed IE lines have been determined to identify 
undulations and peaks of P and T waves. Undulations 
are very small P waves having magnitudes much 
smaller than those of the normal P waves. These are 
also known as the fibrillatory waves [1, 2, 45, 46]. 


In this work, a technique has been developed to detect 
P and T peaks between two consecutive QRS 
complexes based on an algorithm, which follows a set 
of decision logic rules. This technique discriminates 
well between P and T waves so that the next set of 
useful diagnostic features pertaining to the P wave and 
its relationship to the QRS complex is easily 
determined. 


The P wave parameters like P peak location and 
amplitude and P wave onset and offset location and 
amplitude are determined on the basis of slope 
computations performed on the filtered non-QRS 
segments. As per the CSE recommendations, P onset 
has been used as a reference point to determine the 
base level for P wave amplitude measurements. 


Decision Logic Rules to Discriminate Between P 
and T Waves 


It is well known that in arrhythmic signals, P and T 
waves have inconsistent locations w.r.t. the QRS 
complex in successive R-R intervals. For example, 
abnormal rhythms like  bigeminy,  trigeminy, 
ventricular tachycardia or multiple arrhythmia signals 
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display irregularities in locations of P and T waves. In 
such cases, it is very difficult to detect whether a wave 
is a P wave or a T wave in an R-R interval. Absence of 
universally accepted criteria for detecting the 
difference between P and T wave, and the need for 
designing an automatic real time procedure that would 
distinguish between the two waves has motivated 
towards the development of decision logic rules for the 
purpose. 


Different decision logic rules have been developed for 
discriminating between P and T peaks in an R-R 
interval after an exhaustive study of arrhythmia signals 
in the MIT/BIH arrhythmia database and arrhythmic 
ECG signals available in literature [1-5, 24, 46]. 
Consultations with two experienced cardiologists were 
also made in this regard. 


With the help of these rules /conditions, presence, 
polarity, position and number of P/T waves /peaks, 
w.r.t. midpoint in an R-R interval can be detected. All 
these conditions are represented by a parameter termed 
as ‘FLAG’. 


The position of midpoint in an R1-R2 interval is 
computed as equal to (QRS1.¢ + QRS2,,)/2, where 
QRS 1; is offset of the QRS1 complex and QRS2,, is 
onset of the QRS2 complex. A positive peak (P/T) is a 
peak above (upright wave) the detected JE line, and a 
negative peak (P/T) is one below (downright wave) 
this line. 


These decision logic rules are summarized in the look 
up Table 2. The first column in this table is designated 
as FLAG parameter formulated specifically in the 
present study that takes on the values from 1 to 24, 0 
and 100, which govern all the decision rules 
summarized in the remaining columns. 


For instance, when FLAG = 20 the decision logic rule 
is: IF there is one positive peak in the R-R interval 
(RR,), its position w.r.t. the midpoint of the interval is 
towards the right, the consecutive R peaks (R; and R2) 
being +ve in the interval, S peaks (S, and S,) having a 
magnitude greater than 0.24 mV w.r.t. base level, and 
heart rate (HR) greater than 118 bpm; THEN, this 
peak represents a T wave. 


Diagnostic Features Extracted/Deduced from ECG 
Signal 


The 37 useful diagnostic features that have been 
extracted /deduced for arrhythmia detection using the 
above methodology are as follows: 
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Features Extracted from QRS Complex 


R peak sample number, R peak amplitude (mV), R 
peak angle, Q peak sample number, Q peak amplitude 
(mV), S peak sample number, S peak amplitude (mV), 
S peak angle, R’ peak sample number, R’ peak 
amplitude (mV), QRS onset, QRS onset amplitude, 
QRS offset, and QRS offset amplitude are extracted 
from QRS complex 


Features Deduced from QRS Complex 


QRS width, R-R interval, Heart rate (bpm), VAT 
(ventricular activation time), amplitudes of Q, R, R’ 
and S peaks w.r.t. QRS onset, amplitude of QRS offset 
w.r.t. amplitude of QRS onset for a given QRS 
complex, QRS configuration (Table 1) are determined 
from the extracted features of the QRS complex. 


Features Extracted from R-R interval 


Detection of isoelectric line between two consecutive 
QRS complexes in the R-R interval, presence of P/T 
waves, polarity (+/-) of P/T waves, and position of P/T 
waves w.r.t. mid-point of R-R interval are the extracted 
features from R-R interval. All these features are 
represented by a parameter termed as ‘FLAG’, which 
is assigned integer values (+ve) that indicate the 
presence, polarity, position and number of P/T 
waves in an R-R interval w.r.t. isoelectric line 
detected between two consecutive QRS complexes 
and w.r.t. mid-point of R-R interval. 


Features Extracted in R-R interval from P wave 


P peak sample number, P peak amplitude (mV), P 
wave onset, P wave offset, P wave onset amplitude 
(mV), and P wave offset amplitude (mV) are the 
extracted features in R-R interval from P wave 


Features Deduced from R-R intervals 


P-R intervals, P-R segments, P-P intervals, amplitude 
of P peak w.r.t. its onset, amplitude of P wave offset 
w.r.t amplitude of P wave onset for a given P wave, 
and AVCR (Atrio-Ventricular Conduction Ratio) are 
deduced from R-R interval parameters. 


EVALUATION OF DEVELOPED METHOD 
USING SIMULATED SIGNALS 


The algorithms developed for the detection of QRS 
complex and its associated parameters, and for the 
detection of P wave and its associated features, 
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described above, have first been tested using simulated 
ECG signals. Testing any algorithm by applying it to a 
simulated signal helps to create an approximate 
software scheme, which when used for real time data 
analysis aids in the selection of optimised values of 
different parameters related to the logic and sequence 
of data interpretation. 


Keeping this in mind, artificial ECG signals, sampled 
at 360 Hz, have been simulated for a number of 
abnormal beats and rhythms, as per the requirement of 
an arrhythmia type regarding the relative locations and 
amplitudes of the P waves and the QRS complexes. 
Signals that have been simulated conform to the 
morphological variations present in the MLII ECGs of 
the MIT/BIH arrhythmia database. 


While the simulated ECG signals were designed and 
analysed to test the algorithms for feature extraction, a 
number of detailed calculations had to be carried out to 
adjust the tolerance limits and threshold criteria for the 
37 useful diagnostic features listed in Section 3.2. 


After the simulated data test runs were conducted, a set 
of algorithms was obtained which provided a powerful, 
dynamic and automated ECG analysis model that was 
subsequently validated by the real-time ECGs acquired 
from standard data. 


VALIDATION OF DEVELOPED METHOD 
USING STANDARD ECG DATA 


Validation using CSE Data 


If an algorithm detects ECG features within the 
tolerance limits prescribed by the CSE group [39], then 
the influence of errors in measurement of ECG 
parameters by the algorithm does not affect the 
diagnostic significance of these parameters. 


To facilitate the validation of an ECG analysis 
algorithm, 25 records containing 8-10 s of 15 lead 
ECG data sampled at 500 Hz have been provided with 
measurement results and recommended tolerances for 
P onset, P offset, QRS onset, QRS offset and T end in 
CSE Data Set-3 [34,39]. 


For lead II, the recommended tolerance for P onset is 
8.0 ms, P offset is 12.8 ms, QRS onset is 7.8 ms, and 
QRS offset is 12.4 ms. In terms of samples, these 
tolerances are 4.0, 6.4, 3.9 and 6.2 respectively, for 
CSE records sampled at 500 Hz, and, 2.88, 4.608, 
2.808 and 4.464 respectively for the MIT/BIH records 
sampled at 360 Hz. 
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In view of the above, the QRS detection algorithm to 

detect R peak (location) as the QRS fiducial, was 

validated using all the 25 records for lead II of the CSE 

Dataset-3. Lead II ECG signals were specifically 

chosen, since the algorithms in the present study have 

been used to analyse modified lead II ECGs of the 
MIT/BIH data, for arrhythmia analysis. 


The performance of the algorithm for QRS detection in 
terms of R peak fiducial was found to be satisfactory 
for all the 25 records (record nos. 1, 6, 11, 16, 21, 26, 
31, 36, 41, 46, 51, 56, 61, 66, 71, 76, 81, 86, 91, 96, 
101, 106, 111, 116 and 121), and all the QRS 
complexes were detected successfully. 


The composite method of feature extraction was also 
validated for computational accuracy of P onset, P 
offset, QRS onset and QRS offset using lead II ECG 
signals of records 16, 21, 31, 66, 96, 101 and 106 of 
Data Set-3. The numerical values of these four 
parameters, computed by the software developed, were 
found to be within the recommended tolerances for all 
the 85 beats present in these 7 records. 


Using the measurement results of 85 beats present in 
the 7 records, the computed values for R peak 
amplitudes (w.r.t. the baselevel), R-R interval, FLAG 
parameter, P-R interval, P-R segment, P-P interval, and 
AVCR were also validated. 


Validation using MIT/BIH Data 


The process of validation was further carried out for 
the detection of QRS complexes, in terms of R peak 
fiducials from the annotated records of the MIT/BIH 
database. All the records are annotated at R peak 
points, and hence provide true locations of QRS 
complexes. 


The QRS detection algorithm was validated using 
32,800 beats acquired from records 100, 101, 103, 105, 
109, 113, 114, 116, 209, 213, 214, 219, 231, 232, and 
234. Except 14 false detections (9 false positive and 5 
false negative), 32,786 beats were detected correctly, 
wherein 7, 21, 292 QRS features out of 721,600 
extracted QRS features were found to be correct and a 
detection accuracy of 99.96% within the tolerance 
limits recommended by the CSE Working Party. The 
QRS complexes detected as false positive were 
isolated QRS-like artifacts. 


Too much tilt in an IE line is undesirable as it leads to 
erroneous detection of P and T peaks. The amplitude 
band defined for identification of R peak point has a 
special significance in this context (refer to Fig. 4). It 
decides the range in which the R peak amplitudes lie in 
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the ECG signals with a baseline wander from 0.05 to 
0.25 Hz. Due to this band, the IE line computed in an 
R-R interval has a limited tilt that does not effect the 
detection of P and T peaks. Starting from the P peak 
point determined w.r.t. the IE line, P onset and P offset 
are detected using the slope criteria only for the filtered 
signals. 


Therefore, tilted IE lines do not influence their 
computations, and therefore, these P wave boundaries 
are detected within the recommended tolerances given 
by the CSE Working Party. These computations were 
verified manually by observation of the plots of the 
signal under analysis. The actual P peak amplitude is 
computed w.r.t. P onset, which is taken as the base 
level for all P wave measurements. In cases of serious 
baseline wander, to compute P waves and their 
parameters in the time domain, techniques for 
suppression or removal of baseline wander need to be 
employed. 


Authors in [18] have reported that onsets and offsets of 
P and T waves may not have accuracy essential for 
diagnostics when these are influenced by heavy noise 
or baseline wander or when their magnitudes are very 
small. 


The present method is effective for most of the records 
of the MIT/BIH database that are contaminated by 
noise within acceptable levels. Very noisy signals are 
exceptional data segments and in the authors’ 
knowledge have been known to create difficulty for 
other techniques of ECG analysis as well. As also 
pointed out by Gang et al. [20], reliability of any 
technique may not be judged only on the basis of its 
ability to analyse very noisy signals. 


Performance Evaluation of the Composite Method 


The composite method of feature extraction was 
applied to detect R peaks, QRS features like QRS 
onset, QRS offset, QRS width etc., IE line, P/T peaks 
in an R-R interval and P wave features like P onset, P 
offset, P wave magnitude etc. from 8000 beats of 
records 100, 101, 103 and 219. Consistency in the 
values of FLAG parameter was also assessed using 
these beats. 


The software extracted all the 37 features to an 
accuracy that was well within the recommended 
tolerance limits. Verification of the results obtained by 
the method was performed on the basis of detailed 
manual measurements of the ECG plots constituted by 
8000 beats: 2200 beats of record 100, 1800 beats of 
record 101, 2000 beats of record 103, and 2000 beats 
of record 219 of the MIT/BIH database that were 
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successfully identified as QRS complexes by the 
software. 


An unclassifiable beat in a data segment of Record 101 
was detected as QS wave. Varying morphology of the 
QRS complexes in record 219 did not pose any 
problems in the satisfactory detection of QRS features. 
Presence of undulations as fibrillatory waves was also 
identified in the rhythms of this record. 


The abnormal QRS complexes in this record have 
large amplitudes but have slopes that do not fall within 
55% of the highest slope of a rising R wave (slope 
threshold criteria for QRS detection) in the rhythm 
under analysis (2500 sample points are analysed at a 
time). To detect these complexes, a suitable width 
criteria (if the QRS width is less than or equal to 1.5 
times the stipulated maximum width of normal QRS 
complexes then the tentative R peak is marked as a 
valid QRS complex) was employed that also 
discriminate these large complexes from comparable T 
waves. 


After the extensive process of testing and validation of 
the composite method of feature extraction over a 
considerably large amount of standard ECG data the 
software was applied again to extract all the 37 
features from each of the 12 s pathological data 
segments of records 100, 119, 201, 202, 205, 209, 219, 
232 and 233 (approx. 250 beats). These ECG signals 
were not used earlier for analysis. Various 
characteristics found in these data segments are as 
follows: 


Normal Rhythm, record 100; Ventricular bigeminy 
(with uniform PVCs having +ve R wave), Ventricular 
trigeminy (with uniform PVCs having +ve R wave), 
record 119; Pure atrial fibrillation, record 201; Sinus 
bradycardia, record 202; R-on-T premature ventricular 
complex, Slow ventricular tachycardia, record 205; 
Supra-ventricular tachycardia, record 209; Atrial 
fibrillation + premature ventricular complex + T wave, 
record 219; Sinus bradycardia + right bundle branch 
block + A-V block I, record 232; and, Isolated 
premature ventricular complex (QS wave) + atrial 
premature complex, Isolated premature ventricular 
complex with prominent S peak (+ve R wave), record 
233. 


Discussion on Results 


In all the 12 s pathological recordings, a large variation 
is present in the morphology of the QRS complexes. 
Also the P and T waves have different magnitudes and 
polarities in these signals. The signals were analysed 
with noise within acceptable levels and a baseline drift 
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in the range from 0.05 Hz to 0.25 Hz. All the 37 
features were successfully detected in these signals 
with 100% accuracy, i.e. the accuracy of detection was 
satisfactory enough to classify these parameters into 
different diagnostic categories. 


Once the correct value of the FLAG parameter was 
identified by the program, computation of P wave 
features was relatively easy. Width of the QRS 
complex is widely used to identify premature 
ventricular complexes, and its correct determination 
has also been of much significance. 


The maximum value of QRS width for practically all 
normal QRS complexes detected in the signals was 
found to be 0.115 s, although its upper limit has been 
specified as 0.12 s in general, and 0.10 s, 0.11 s in 
particular [1-3]. 


Some wide and abnormal QRS complexes were found 
to have their widths as low as 0.117 s. In record 202, 
the width of few normal QRS complexes was found to 
be 0.122 s. However, to detect bradycardia in record 
202, the QRS width criteria is not required (heart rate 
computations are sufficient), and in general the 
maximum limit of normal QRS width was fixed at 
0.115 s. 


Since the signals subjected to analysis have been 
shifted upwards by the baseline (0 mV), to obtain only 
positive amplitudes of ECG waves, all the amplitudes 
have been initially computed from this shifted 
baseline. Subsequently, the exact amplitudes of QRS 
wave were computed w.r.t. to the QRS onset as the 
base level while those of P wave were computed w.r.t. 
P onset as the base level. 


More than 8250 R-R intervals were filtered and 
analysed in a number of normal and abnormal rhythms 
of the MIT/BIH database. During this process of 
validation, the extracted/deduced features were 
measured in the filtered signals and manually verified 
against those in the unfiltered ECGs. The results of 
two representative abnormal rhythms for MIT/BIH 
Rec. No. 205 with signal duration 1677.989 — 1684.3 s, 
and, MIT/BIH Rec. No. 202 with signal duration 
1107.333 — 1113.997 s, are shown in Figs. 8-9, Tables 
3-4 and Figs. 10-11, Tables 5-6, respectively. 


The length of the data segments shown has been 
restricted due to limitation of space. Figures 8-9 and 
10-11 show the unfiltered signal and the signal with 
filtered non-QRS portions for Rec. Nos. 205 and 202, 
respectively. Tables 3-6 present the entire features 
extracted /deduced i.e., features pertaining to QRS 
wave and R-R interval (refer Sec. IIIB), for the above 
two records (for symbols used in Tables 3-6, refer 
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Table 7). For example, Table 3 gives the 22 features 

extracted/deduced pertaining to each of the 9 QRS 

waves present in the signal duration 1677.989 — 1684.3 

s of Rec. No. 205, while Table 4 gives the 15 features 

extracted/deduced pertaining to each of the 8 R-R 
intervals present in the above record. 


The automated extraction of diagnostic features, using 
the composite methodology, can be employed by a 
suitable classifier such as a knowledge-based expert 
system for easy and reliable interpretation of single as 
well as multiple arrhythmic events in an ECG signal. 
For example, using the above extracted features in the 
knowledge-base of an expert system, it can be reliably 
interpreted that in Rec. No. 205, the 7 QRS wave 
shows isolated R-on-T Premature Ventricular Complex 
(PVC) with downright R peak, and the interval 
between 7" and 8" R peaks shows Full Compensatory 
Pause (FCP). 


ADVANTAGEOUS FEATURES OF THE 
DEVELOPED METHODOLOGY 


A number of computer-based techniques are reported 
in literature for feature extraction of ECG signals. As 
the QRS complex is more distinguishable and easier to 
detect in comparison with other ECG waves, the R- 
peak point is generally used as a reference to reliably 
detect the QRS wave; on the basis of QRS detection, 
other prominent features are extracted from the ECG 
signal [6, 47]. 


Various methods developed so far for detection of 
QRS complex can be broadly classified into four 
categories, viz., syntactic methods, non-syntactic 
methods, transformative methods and hybrid methods 
[26, 36, 47-49]. 


The diversity of algorithms developed for QRS 
detection and extraction of the P wave, is too large to 
compare with the algorithms developed in the present 
study for detection of the QRS and P waves. Therefore, 
the comparison has mainly been restricted to 
distinctive QRS complex and P wave detection 
methods, applied to signals of the MIT/BIH database. 


Pan and Tompkins [13] developed a real-time QRS 
detection algorithm to reliably recognize QRS 
complexes based upon analyses of slope, amplitude, 
and width. Hamilton and Tompkins [14] implemented 
another real-time QRS detection algorithm using a 
preprocessor stage to produce event vectors consisting 
of peak heights and time of peak occurrence, on the 
basis of the previous method [13], from the processed 
ECG signal. 
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In the next stage of ECG processing, complex decision 
tules were applied progressively to these event vectors 
to evaluate the quantitative results of various decision 
tule components. In both the methods, P wave 
detection has not been addressed. 


Similarly, Lin and Chang [7] used linear prediction for 
QRS feature extraction only. Xue et al. [15] developed 
a multilayered neural network based adaptive matched 
filter to detect QRS complexes from ECG signals. The 
technique does not report the diagnostic value of 
detected QRS complexes’ regarding different 
arrhythmias in context to differing morphologies. 


Murthy and Prasad [16] developed a pole-zero models 
to measure P, QRS and T waves. However, no method 
has been given to distinguish a P wave from a T wave 
in real-time ECG signals. 


Woolfson [8] applied Kalman filtering to short 
stretches of ECG data for computing the time varying 
spectra of the R-R intervals. Variations in the 
morphology of the QRS complex and detection of the 
P wave, two important parameters necessarily needed 
to interpret arrhythmias, have not been addressed in 
this study. 


Poli et al. [17] used genetic algorithms to detect QRS 
complexes in the ECG signals of the MIT/BIH 
database. However, the diagnostic importance of the 
detected QRS complexes has not been mentioned in 
the study. Li et al. [18] used an algorithm based on the 
wavelet transform to first detect the QRS complex, 
then the T wave, and finally the P wave from single 
channel records. 


As, acknowledged by the authors, onsets and offsets of 
P and T waves were not detected to an accuracy 
required for reliable arrhythmia interpretation. 
Sahambi et al. [12] used wavelet transform to obtain 
multiscale analysis for timing characterization of the 
ECG, for detecting the QRS complex, and, P and T 
waves with positive and negative polarity. The study 
has presented the results only for record 106 of the 
database. 


Wieben et al. [9] used the concept of filter bank to 
implement time-frequency analysis of MLIT ECG 
records. The authors also analysed these ECG signals 
to obtain heart rate and QRS morphological features. 
The filter bank features, and the heart rate and QRS 
features were used to classify only premature 
ventricular complexes. 


Gang et al. [20] used linear-approximation distance- 
thresholding compression with a _ backpropagation 
neural network to detect normal and abnormal QRS 
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complexes from eight records of the MIT/BIH 

database. However, the method has not been used to 

detect P and T waves in these ECG signals. Bartolo et 

al. [10] have developed QRS detection software based 
on QRS template matching process. 


Record segments with poor signal quality, where most 
of the QRS complexes were visually indistinguishable 
from background activity, were excluded from 
analysis. As suggested by the authors, further 
improvement in discriminating normal and abnormal 
premature beats required detection and analysis of the 
P wave that was not attempted in the study. 


Osowski and Linh [11] statistically transformed QRS 
complex to obtain a suitable feature vector. The feature 
extraction method does not detect the P wave that is of 
utmost importance to classify arrhythmias. Botter et al. 
[19] used a_ backward-error-propagation neural 
network with asymmetric basis functions to extract the 
features of P waves as a set of nine parameters. 


The method does not address detection of the QRS 
complex, in the absence of which arrhythmias can 
hardly be classified, except the atrial fibrillation and 
the atrial flutter rhythms. Owis et al. [21] presented a 
study of the non-linear dynamics of ECG signals for 
arrhythmia characterization. The results indicated that 
the discrimination between different arrhythmia types 
is difficult using features obtained by this study. 


Shyu et al. [22] used quadratic spline wavelets to 
detect QRS complexes from seven records of the 
MIT/BIH database. P wave detection has not been 
addressed in the study. Chazal et al. [23] presented 
automatic classification of heart beats using ECG 
morphology and heart beat interval features. In this 
work, QRS _ fiducials provided by the standard 
annotations have been directly utilized to obtain 
different ECG features. 


R peak detection has not been performed. Illanes- 
Manriquez and Zhang [49] detected QRS onset and 
offset in ECG signals by computing an indicator 
related to the area covered by the QRS complex 
envelope. However, P wave detection was not 
addressed in the study. 


On comparing the composite method, developed in the 
present study, with the methods used by other 
researchers discussed above to detect QRS and P 
waves from ECG signals of the MIT/BIH database, 
four unique contributions of the developed composite 
method emerge: 


e ~=Extraction of QRS complexes in such a way that 
their original morphology is preserved for accurate 
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determination of several features of diagnostic 
importance, and detection of P waves in a manner 
that these can be well discriminated from T waves 
in different abnormal rhythms; 


Demonstration of the inherent strength of the 
method for signal classification that is not so 
transparent in other known techniques applied to 
the analysis of a variety of real-time normal and 
abnormal beats and rhythms; providing an 
amenable technique that can be appreciated and 
modified by a cardiologist; 
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A simple and understandable technique (devoid of 
any ‘black-box’ effect) that is more practical in a 
clinical environment, and provides a powerful 
medical information processing system that can 
be successfully employed to train medical 
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Fig. 1. Normal ECG cardiac cycle 
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Fig. 2. Flow diagram depicting the steps followed in the composite method of feature extraction 
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Fig. 3. Slurs (A, C, E, G) and notches (B, D, F, H) in component waves (Q-R-S-R’) of the QRS 
complex 
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Fig. 4. Flowchart for recognition of QRS wave and characteristic point identification 
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Fig. 5. Amplitude responses of 6th order low pass Butterworth Filter with cut off frequency of 10 Hz 
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Fig. 6. Flow chart for detection of IE line 
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Fig. 7. Flow chart for recognition of P wave and characteristic point identification 
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Fig. 9. MLII signal with filtered non-QRS portions (MIT/BIH Rec. No. 205) 
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Fig. 8 MLII signal with filtered non-QRS portions (MIT/BIH Rec. 202) 
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Table 1 Different configurations of the QRS complex (MLII) 


QRS morphology Configuration designation 


R 
| R’ QRSR’ 
Qs 
R 
QRS 
Qs 
R 
| 7 
Q 
R 
| \* RSR’ 
S 
R 
| RS 
S 
R 
AL 
VO . 


Single negative peak 
(R wave absent) 
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Table 2. Look-up table of decision logic rules developed for discriminating between P and T peaks/waves in an R-R interval 


No. of peaks present in R-R 


FLAG ai Polarity (+ve/-ve) of peaks in R-R interval Position of peaks with respect to midpoint in R-R interval (L, M, Additional criteria to determine peak/wave type (P or T) in R-R interval Discrimination of peak/wave type (P or T) present in R-R 
(RR,)* R)* (RR,)* interval 
2"7 peak 2"™ peak 
1 +ve +ve - L R - HR = 60-100 bpm in RR, T P - 
2 1 +ve - - L/M - - HR>118 bpm in RR,; RR, with +ve R; & +ve R2 both with S waves Tr - - 
3 2 -ve +ve - L L - - P £ - 
4 2 +ve +ve - L L - - T P - 
5 2 -ve +ve - L R - HR>118 bpm in RR,;; RR; with +ve R; & +ve R2 both with S waves P Tr - 
6 1 +ve - - R/M - - HR = 60-100 bpm in RR,; RR; with +ve R; & +ve Ro P - - 
7 1 -ve - L - - RR, with +ve R;, Ro T with undulations after it 
8 1 +ve - - L/M RR, with -ve R; & -ve R; & no S waves T - - 
9 2 -ve +ve LL. R - T P : 
10 3 -ve +ve +ve L L R RR, with +ve Ri, R2; S;, S2 > 0.24 mV P T P 
11 1 +ve/-ve - - M/R - - RR,<(RR,1 or RRq1) Tr - - 
12 1 -ve - - R - - Peak amplitude >= half the amplitude of higher R peak in RR, T - - 
13 1 +ve R HR = 60-100 bpm in RR; P with undulations before it 
14 3 -ve +ve +ve 1 L/M R HR = 60-100 bpm in RR, Di-phasic T P 
15 2 +ve/-ve +ve/-ve L/M/R L/M/R Peaks are extremums of sawtooth waves P/flutter wave P/flutter wave - 
16 3 or more +ve +ve +ve L/M/R L/M/R L/M/R HR < 120 bpm in RR;; S;, S2 < 0.24 mVor absent in RR, All the peaks are P peaks 


Int L: left of mid-point; M: at mid-point; R: right of mid-point of R-R interval; QRS(W): width of QRS complex; RR;,: interval (in sec) between two consecutive R peaks, say R,; & Ry; RR,.;: R-R interval preceding RR,; RR,,,: R-R interval succeeding RR, ; 
HR : heart rate in bpm; ‘-’ symbolizes absence of respective peak/ wave in corresponding R-R interval. 


(contd...) 
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No. of peaks present in R- Polarity (+ve/-ve) of peaks in R-R Position of peaks with respect to midpoint in R-R | Additional criteria to determine peak/wave type (P or T) Discrimination of peak/wave type (P or T) 
FLAG R interval interval (RR,)* interval (L, M, R)* in R-R interval (RR,)* present in R-R interval 
1“ peak 2" peak 3" peak 1* peak 2™ peak 3" peak 1“ peak 2" peak 3" peak 


RR, with —ve Ro, or +ve R2 with prominent S wave; 


RR,<(RR:1 or RR); ~HR>100 bpm in RR; 


RR, with +ve R; & —ve Ro; peak within 0.15 s of R2 


RR, = RRi1 


RR, with +ve Rj, Ro; S;, S2 > 0.24 mV; HR>100 bpm in RR, 


RR, with +ve Ry, Ro; Si, S2 > 0.24 mV; HR>100 bpm in RR, False +ve 


RR, with +ve Rj, R2; QRS1(W) & QRS2(W) > 0.115 s 


RR, with +ve Ry, Ro; S1, S2 < 0.24 mVor absent; QRS1(W) & P with undulations after it 
QRS2(W) < 0.115 s 


RR, with +ve R; & -ve Ro; peak not within 0.15 s of Ro 


RR, with +ve R; & -ve Ro No peaks present 


None of the above conditions fulfilled 


International Journal of BioSciences and Technology (2009), Volume 2, Issue 7, Page(s):96-125 
118 


Table 3. QRS features computed from MIT/BIH rec. no. 205 
Duration: 1677.989-1684.3 sec 


QRS Features 


1678.264 
1.200 


1678.292 
2.495 
1.295 
2.446 

1678.311 
1.135 
0.065 
1.962 

1678.328 
1.405 
0.205 

1678.353 
1.285 
0.085 
0.089 
0.028 
RSR' 


Samp (mV) 
Sang (deg) 


QRS config 


1679.006 
1.200 


1679.031 
2.510 
1.310 
2.412 

1679.053 
1.135 
0.065 
2.247 

1679.069 
1.370 
0.170 

1679.103 
1.280 
0.080 
0.097 
0.025 
RSR' 


1679.703 
1.250 
1679.719 
1.155 
0.095 
1679.750 
2.520 
1.270 
2.343 
1679.767 
1.115 
0.135 0.130 
1.935 2.407 

- 1680.528 
- 1.395 
. 0.110 
1679.781 1680.550 
1.335 1.290 
0.085 0.005 
0.078 0.108 
0.047 0.050 
QRS QRSR' 


1680.442 
1.285 
1680.464 
1.220 
0.065 
1680.492 
2.510 
1.225 
2.315 
1680.508 
1.155 


1681.172 
1.280 
1681.192 
1.200 
0.080 
1681.219 
2.470 
1.190 
2.505 
1681.239 
1.150 
0.130 
2.609 


1681.258 
1.375 
0.095 
0.086 
0.047 
QRS 


1681.939 1682.311 
1.180 1.525 


1681.967 1682.347 
2.460 0.000 
1.280 1.525 
2.460 6.173 

1681.986 
1.095 
0.085 
1.960 

1682.003 
1.365 
0.185 

1682.022 
1.250 
0.070 
0.083 
0.028 
RSR' 


1682.447 
1.500 


* Refer list of symbols (Table 7) for description of variables listed in QRS features. 
symbolizes absence of respective QRS feature in the corresponding QRS complex. 


Table 4. R-R interval features computed from MIT/BIH rec. no. 205 
Duration: 1677.989-1684.3 sec 


R-R Interval R1-R2 
Features 


0.739 
81.191 
1 
1678.853 
1.272 
1678.900 
1.323 
0.051 
1678.956 
1.243 
0.029 
0.153 
0.050 
0.717 
11 


Pamp (mV) 
Pf (s) 

Pfa (mV) 
Pfamp (mV) 
P-R int (s) 
P-R seg (s) 
P-P int (s) 
AVCR 


R2-R3 


0.719 
83.449 
1 
1679.569 
1.265 
1679.617 
1.316 
0.051 
1679.675 
1.227 
0.038 
0.134 
0.028 
0.741 
111 


R3-R4 


0.742 
80.863 
1 1 
1680.303 
1.286 
1680.358 
1.353 
0.067 
1680.417 
1.271 
0.015 
0.139 
0.025 
0.731 
1:11 


0.727 


1.297 


1.343 
0.046 


1.261 
0.036 
0.130 
0.025 
0.747 
11 


R4-R5 


82.531 


1681.042 


1681.089 


1681.147 


R5-R6 


0.748 
80.214 


{ 


1:1 


* Refer list of symbols (Table 7) for description of variables listed in R-R interval features. ‘-’ 


1681.783 
1.264 
1681.836 
1.320 
0.056 
1681.897 
1.231 
0.033 
0.156 
0.042 


R6-R7 


0.380 
157.895 
0 


symbolizes absence of respective feature in the corresponding R-R interval. 
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1683.400 


1684.136 
1.180 


1683.414 - 


1683.442 


1683.461 


1684.164 
2.520 
1.340 
2.359 

1684.183 
1.105 
0.075 
1.682 


1683.475 


1683.492 


R7-R8 


1.095 
54.795 


1 


1683.269 
1.236 
1683.317 
1.281 
0.045 
1683.367 
1.223 
0.013 
0.131 
0.033 
0.716 


111 


1684.197 
1.375 
0.195 
0.061 
0.028 

RS 


R8-R9 


0.722 
83.102 
1 
1683.975 
1.272 
1684.033 
1.332 
0.060 
1684.094 
1.249 
0.023 
0.161 
0.042 
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Table 5. QRS features computed from MIT/BIH rec. no. 202 
Duration: 1107.333-1113.997 sec 


1107.725 
0.130 


1108.889 
0.120 


1110.028 
0.140 


1111.139 = 1112.325 
0.225 0.210 


1107.775 1108.942 1110.081 1111.186 1112.372 
2.165 2.070 2.100 2.200 2.140 
2.035 1.950 1.960 1.975 1.930 
Rang (deg) 3.577 3.256 3.052 2.608 2.783 
Ss (s) 
Sa (mV) 


Sang (deg) 

R's (s) 

R'a (mV) 
R'amp (mV) 
QRSf (s) 
QRSfa (mV) 
QRSfamp (mV) 


1107.833 1108.981 1110.114 1111.214 1112.403 
0.060 0.175 0.215 0.245 0.255 
0.070 0.055 0.075 0.020 0.045 

QRS(W) (s) 0.108 0.092 0.086 0.075 0.078 

VAT (s) 0.050 0.053 0.053 0.047 0.047 

QRS config R R R R R 


* Refer list of symbols (Table 7) for description of variables listed in QRS features. 
symbolizes absence of respective QRS feature in the corresponding QRS complex. 


Table 6. R-R interval features computed from MIT/BIH rec. no. 202 


Duration: 1107.333-1113.997 sec 


R-R Interval R1-R2 
Features 


1.167 
51.414 
1 
1108.675 
0.071 
1108.758 
0.172 
0.101 
1108.825 
0.069 
0.002 
0.214 
0.064 
1.150 
1:1 


* Refer list of symbols (Table 7) for description of variables listed in R-R interval features. ‘- 


R2-R3 


1.139 
52.678 
1 
1109.836 
0.090 
1109.908 
0.177 
0.087 
1109.969 
0.099 
0.009 
0.192 
0.059 
1.100 
1:1 


R3-R4 


1.105 
54.299 
1 
1110.947 
0.199 
1111.008 
0.274 
0.075 
1111.078 
0.152 
0.047 
0.192 
0.061 
1.175 
1:1 


R4-R5 


1.186 
50.590 
1 


1112.111 
0.158 
1112.183 
0.246 
0.088 
1112.256 
0.141 
0.017 
0.214 
0.069 
1.206 
1:1 


symbolizes absence of respective feature in the corresponding R-R interval. 
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1113.525 
0.150 


1113.578 
2.245 
2.095 
2.836 


1113.6114 
0.195 
0.045 
0.086 
0.053 

R 


1.206 
49.751 
1 
1113.308 
0.133 
1113.389 
0.218 
0.085 
1113.456 
0.119 
0.014 
0.217 
0.069 
1:1 


‘° 


WEST 


e The technique extracts the different characteristic 
points of the QRS complex and the P wave (QRS- 
onset, QRS-offset, P-onset, P-offset) within the 
tolerance limits recommended by the CSE 
Working Party. The composite method has the 
potential to detect 37 diagnostic features used by 
cardiologists to classify arrhythmias in lead II 
ECG signals, provided these signals have 
distinguishable QRS complexes and P waves 
contaminated by any type and amount of 
background noise, but with a baseline wander not 
exceeding 0.25 Hz. 


CONCLUSIONS 


ECG signals exhibiting irregular waveforms are 
indicators of arrhythmic events that lead to life 
threatening rhythms, such as ventricular fibrillation, 
ventricular tachycardia and atrial fibrillation. For 
example, R-on-T PVC generally leads to serious 
ventricular tachycardia resulting in deadly ventricular 
flutter /fibrillation. Detecting a single arrhythmia is 
much simpler than detecting a combination of 
arrhythmias (in the same ECG signal), as a small 
number of features extracted is sufficient for 
interpretation. In real situations, as observed in the 
standard ECGs, abnormal rhythms and beats exist in 
combinations of two or more. This overlapping feature 
characteristic of arrhythmias makes it challenging to 
obtain adequate feature detection and classification. 


In this paper, a composite methodology has been 
developed to perform cycle-by-cycle feature extraction 
from the given signal to enable computer-aided 
interpretation of different arrhythmias. This helps the 
cardiologist while diagnosing one arrhythmia at a time, 
or a combination of two or more arrhythmias 
coexisting in the signal. The developed technique has 
been applied to different signals of MIT/BIH Database. 
The QRS detection algorithm was validated using 
32,800 beats of several records of the MIT/BIH 
database and a detection accuracy of 99.96% was 
achieved within the tolerance limits recommended by 
the CSE Working Party. The composite method for 
detection of both P and QRS wave features has been 
validated using all the 25 records for lead II of the CSE 
Dataset-3, before being applied to approximately 8000 
beats of the MIT/BIH database, wherein nearly 
2,96,000 P wave and QRS complex features were 
extracted and manually verified. 
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The technique eliminates most of the shortcomings of 
the existing methods for ECG feature extraction, and 
leads to the detection of a number of useful diagnostic 
features related to the depolarisation sequence of the 
heart. On the basis of this strong set of diagnostic 
features, not only the isolated abnormal beats or 
rhythms can be detected, but also combinations of 
normal /abnormal beats and rhythms, as they appear in 
an ECG signal, can be successfully interpreted. The 
methodology is dynamic and amenable to further 
refinements if needed by a cardiologist to extend its 
potential to detect as many arrhythmias as possible. 
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Table 7. List of symbols 


Symbol Description 


AVCR Atrio-ventricular conduction ratio 


FLAG Integer number (+ve) indicating the presence, polarity, position and number of P/T waves in an interval w.r.t. 
isoelectric line detected between successive QRS complexes and w.r.t. mid-point of RR interval 


HR Heart rate in beats per minute (bpm) 

Pa Amplitude (in mV) of P peak 

Pamp |Pna - Pal 

Pf Location of P wave offset in seconds corresponding to a sample number 

Pfa Amplitude (in mV) of P wave offset 

Pfamp |Pna - Pfa| 

Pn Location of P wave onset in seconds corresponding to a sample number 

Pna Amplitude (in mV) of P wave onset 

Ps Location of P peak in seconds corresponding to a sample number 

P-R int Interval (in sec) between onset of P wave and onset of the following QRS complex 
P-R seg Interval (in sec) between offset of P wave and onset of the following QRS complex 
P-P int Interval (in sec) between onset of P wave and onset of the next P wave 

Qa Amplitude (in mV) of Q peak 

Qamp |QRSna - Qa| 

Qs Location of Q peak in seconds corresponding to a sample number 

QRSf Location of QRS complex offset in seconds corresponding to a sample number 
QRSfa Amplitude (in mV) of QRS complex offset 
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QRSfamp = |QRSna - QRSfa| 


QRSn Location of QRS complex onset in seconds corresponding to a sample number 
QRSna Amplitude (in mV) of QRS complex onset 

QRS Configuration of QRS complex (i.e. QRSR’, QRS, QR, RSR’, RS, R, QS) 
config 


QRS(W) Time duration or width (in sec) of QRS complex 


Ra Amplitude (in mV) of R peak 

Ramp |QRSna - Ral 

Rang Angle of R peak (in degrees) 

Rs Location of R peak in seconds corresponding to a sample number 
R’a Amplitude (in mV) of R’ peak 

R’amp |QRSna — R’a| 

R’s Location of R’ peak in seconds corresponding to a sample number 
R-R int Interval (in sec) between two consecutive R peaks 

Sa Amplitude (in mV) of S peak 

Samp |QRSna - Sal 

Sang Angle of S peak (in degrees) 

Ss Location of S peak in seconds corresponding to a sample number 
VAT Ventricular activation time 
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